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We report the time dependent response of electrical resistivity in the non-magnetic perovskite 
oxide NdNiC>3 in its phase separated state and provide a physical explanation of the observations. 
We also model the system and do an accurate Monte Carlo simulation of the observed behavior. 
While cooling a phase separation takes place in this system below its metal-insulator transition tem- 
perature and in this state the material exhibits various dynamical phenomena such as relaxation of 
resistivity, dependence of resistivity on cooling rate and rejuvenation of the material after ageing. 
These phenomena signal that the phase separated state of NdNiC>3 is not in thermodynamic equilib- 
rium and we conjecture that it consists of supercooled paramagnetic metallic and antiferromagnetic 
insulating phases. The supercooled phases are metastable and they switch over to the insulating 
equilibrium state stochastically and this can account for the slow dynamics observed in our system. 
We also verify the predictive power of our model by simulating the result of a new experiment and 
confirming it by actual measurements. 
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I. INTRODUCTION 

Usually hard condensed matter is associated with fast 
dynamics, but sometimes we also see slow dynamics, par- 
ticularly in cases such as glasses, spin glasses, phase sep- 
arated systems etc. Dynamics of glasses and spin glasses 
have been intensively studied from the middle of last cen- 
tury, but such studies on phase separated systems started 
attracting a good deal of interest from physicists only 
towards the end of last century. Phase separation al- 
ways seems to be associated with a broadened first order 
phase transition and it has been observed in such sys- 
tems as transition metal oxides, a famous example being 
the colossal magnetoresistance (CMR) manganites, and 
alloys which show martensitic phase transitions.— &&&& 

There has been a large number of studies on the 
phase separated states associated with a broad first or- 
der metal-insulator (M-I) transition in transition metal 
oxides such as VO2, CMR manganites, the rare earth 
nickelate PrNi03, the cobaltite Lai-^ Sr^CoOa etc. They 
exhibit unusual phenomena such as very slow relaxation 
of resistivity& 7 -'& & 10 ' 11 and magnetization , 12 ' 13 rejuvena- 
tion of resistivity after ageing and dependence of physical 
properties on cooling rate^ i 14 ' 15 ' 16 A significant amount 
of work has gone in to study these systems in detail but 
they are not, as yet, well understood even though there is 
a recognition that the phenomena described above orig- 
inate from phase separation. The presence of the slow 
glass-like relaxations in the above systems has led many 
people to consider the phase separated systems as spin 
glasses, cluster glasses and it has been argued that the 
intercluster magnetic interactions are responsible for the 
observed glass-like dynamics in these systems.—^ The 
above picture has been proposed for magnetic systems, 
but similar phenomena have been seen in non-magnetic 



systems as well^ This raises the possibility that a spin 
glass like model may not be the most appropriate to de- 
scribe the dynamics of phase separation in general. There 
has been other attempts to describe these phenomena in 
CMR manganites in terms of a phenomenological model 
where it has been proposed that the phase boundaries 
relax through a hierarchy of energy barriers.^ 2 - Though 
this is a commendable attempt we feel that this model 
does not really capture the essence of the phenomena be- 
cause it does not consider what happens to the system 
while heating. There has been a more recent approach 
to understand the dynamics of phase separated systems 
based on the concept of kinetic arrest of a supercooled 
system into a magnetic glass and its subsequent de-arrest 
on warming i 11 1 15 But there are systems which are non- 
magnetic and non-glassy and it would be very interesting 
to make an attempt to find out the physics behind the 
slow dynamics in such systems because what it might 
reveal may turn out to have wider implications^ ' 17 ' 18 

The perovskite oxide NdNiC>3, which undergoes a tem- 
perature driven M-I transition and has an associated 
phase separated stated we believe, can be taken as a 
model system to probe the dynamics of coexisting phases 
in non-magnetic, non-glassy systems. NdNiC>3 is a clean 
system which does not require any doping to see the 
phase separated state and it is essentially unaffected by 
the application of a magnetic field and it does not show 
any evidence of glass transition . 20 ' 21 Moreover the phase 
separated state exists over a relatively large temperature 
range which makes it easy to study slow dynamics and 
phase separation. 

NdNiC>3 is a member of the series of compounds 
known as rare earth nickelates (RMO3), which are one 
of the few families of perovskite oxides that undergo a 
first order metal-insulator phase transition* 2 ^ These com- 
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pounds crystallize in an orthorhombically distorted per- 
ovskite structure with space group Pbnm- The ground 
state of the nickelates (R^La) is insulating , 22 ' 23 charge 
ordere d 24 i 25 i 26 i 27 i 28 and antiferromagnetic?2& On increas- 
ing temperature from absolute zero these compounds un- 
dergo a temperature driven antiferromagnetic to param- 
agnetic transition, and an insulator to metal transition. 
Below the M-I transition temperature (Tmi), transport 
properties of NdNiC>3 exhibit a large hysteresis and in 
this region, the physical state of the system is phase sep- 
arated and it has been shown to consist of paramagnetic 
metallic and antiferromagnetic insulating phases i 19 i 22 

In this work we use careful time and temperature de- 
pendent resistivity measurements on NdNiC>3 to gain an 
understanding of its phase separated state and its dynam- 
ics. The phase separated state comes into existence on 
cooling the material below its M-I transition temperature 
and vanishes at sufficiently low temperature and does not 
form at all during subsequent heating. We have provided 
a credible physical explanation along with a Monte Carlo 
simulation to describe all the behavior exhibited by our 
system. The simulation results are found to reproduce 
the experimental results quite well. Our results indicate 
that in the phase separated state of NdNi03, the metallic 
phase is present in its supercooled state. The supercooled 
metallic phase is metastable and it relaxes to the insulat- 
ing equilibrium state giving rise to the various observed 
phenomena. 



II. EXPERIMENTAL DETAILS 

Polycrystalline NdNi03 samples in the form of 6 mm 
diameter and 1 mm thick pellets were prepared and 
characterized as described elsewhere.- 3 ^ The preparation 
method uses a high temperature of 1000°C and a high 
oxygen pressure of 200 bar. 

All the temperature and time dependent measurements 
were done in a home made cryostat. To avoid ther- 
mal gradients in the sample during measurement it was 
mounted inside a thick-walled copper enclosure so that 
during the measurement the sample temperature would 
be uniform. It was found that mounting the sample in 
this fashion improved the reproducibility of the time de- 
pendence measurements significantly. A Lakeshore Cry- 
otronics temperature controller was used to control the 
temperature and the temperature stability was found to 
be better than 3mK during constant temperature mea- 
surements. 

Below T M i (~ 200 K) NdNi0 3 is not in thermody- 
namic equilibrium and slowly relaxes because of which 
the experimental data that we get depend on the proce- 
dure used for the measurement. The procedure we used 
was as follows. While cooling we start from 300 K, and 
then record the data in steps of 1 K interval after allowing 
the temperature to stabilize at each point. In between 
two temperature points the sample was cooled at a fixed 
cooling rate of 2K/min. After the cooling run is over 



we wait for one hour at 82 K and then the heating data 
was collected at every one degree interval. The heating 
rate between temperature points was the same as the 
cooling rate used earlier. This cycle of measurements 
was repeated with a different cooling and heating rate of 
0.2 K/min also. 

It was observed that the resistivity above 200 K and 
below 115 K does not show any time dependence, and 
it is also independent of measurement history. Thus to 
avoid the effect of any previous measurements, all time 
dependent experiments in the cooling run were done as 
follows: first take the sample to 220 K (sufficiently above 
200 K), wait for half an hour, then cool at 2.0 K/min to 
the temperature of interest and once the temperature has 
stabilized record the resistance as a function of time. In 
the heating run the time dependent resistivity was done 
in a similar fashion: first take the sample to 220 K, wait 
for half an hour, then cool at 2.0 K/min to 85 K, wait for 
one hour, and then heat at 2.0 K/min to the temperature 
of interest and once the temperature has stabilized record 
the resistance as a function of time. 

The four probe van der Pauw method was used to mea- 
sure the resistivity and standard precautions, such as cur- 
rent reversal to take care of stray emfs, were taken during 
the measurement. We also took care to ensure that the 
measuring current was not heating up the sample. 



III. RESULTS 

Figure Q] shows the electrical resistivity of NdNi03 as a 
function of temperature. The resistivity is multiple val- 
ued, the cooling and heating data differing significantly 
from each other and forming a large hysteresis loop. The 
resistivity plot indicates that NdNiOs undergoes a rela- 
tively sharp M-I transition at about 200 K while heating 
with a width of about 10 K. In contrast, while cooling, 
the resistivity shows a rather broad M-I transition cen- 
tered around 140 K with a spread of about 40 K. Below 
115 K or so, the heating and cooling data merge and the 
logp vs 1/T plot is linear. This indicates that the sample 
is insulating at low temperatures and, if the band gap is 
A, the resistivity should follow the relation 

p(T) = PQ e^ (1) 

Below 115 K both the heating and the cooling data fit 
quite well to this model, with a coefficient of determina- 
tion, R 2 = 0.99955. po and A for the insulating region 
turn out to be 99 mfl cm and 42 meV respectively, which 
are in reasonable agreement with the values previously 
reported^ 

We collected more hysteresis data with different min- 
imum temperatures such as 140 K, 146 K and 160 K. In 
these measurements we cool the sample from 220 K to 
one of the minimum temperatures mentioned above and 
then heat it back to 220 K, both operations being car- 
ried out at a fixed rate of 2 K/min. Loops formed in this 
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fashion are called minor loops and these are indicated by 
the labels a, b and c in Figure [U In the cooling cycle 
all the three minor loops coincide with the cooling curve 
of the full hysteresis loop. In the heating cycle, for loops 
a and b with lower minimum temperatures, the resistiv- 
ity decreases with increasing temperature and joins the 
full loop at 200 K. In the case of loop c, as we increase 
the temperature, the resistivity increases somewhat till 
about 188 K and then it falls and joins the full loop at 
200 K. 

The resistivity also shows a noticeable dependence on 
the rate of temperature change in the cooling cycle as 
shown in the inset of Figure [U The data for the lowest 
curve was collected at 2 K /min and for the middle curve 
at 0.2K/min. The uppermost curve is an estimate ob- 
tained by extrapolating the time dependence data shown 
in Figure [2] to infinite time. We did not see any de- 
pendence on rate of change of temperature when heating 
from 80 K to 220 K. This means that while cooling, below 
the M-I transition temperature (200 K), the system is not 
in equilibrium and on the other hand while heating, the 
system is either in, or very close to, equilibrium. These 
observations are corroborated by the data shown in the 
next Figure. 

A subset of the time dependent resistivity data taken 



2.5 



1.0 





145 K 




147.5 K 








150 K 

^ == ^Z^^~- T4cTk 

. ■ — 1bfaK 


Heating run, 145 K 



1000 2000 3000 

time (s) 



Figure 2: (Color online) Time dependence of resistivity while 
cooling, at various temperatures in the range 140 to 155 K, 
for a period of one hour. The maximum time dependence is 
seen at 145 K which is about 200 %. The curve at the bottom, 
which looks like a horizontal straight line, shows the increase 
in resistivity in a heating run taken at 145 K and the total 
change for one hour in this case is less than 0.2 %. Not all the 
data are shown here to avoid clutter. 
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Figure 1: (Color online) p vs. l/T plot for NdNi0 3 . The 
blue circles represent cooling data and the red squares stand 
for heating data (cooling/heating rate 2.0 K/min). The solid 
line is a least square fit to the band gap insulator model below 
115 K. The curves labeled a, b & c were taken as described in 
the text. The accuracy of the data points is better than 1% 
everywhere. The error does not exceed 0.8 mfi cm anywhere. 
The inset shows p vs l/T for three different cooling rates: 
lower curve (blue stars): 2 K/min, middle curve (red circles): 
0.2 K/min, upper curve (black pluses): infinitely slowly (ex- 
plained in text). The connecting lines are to guide the eyes. 



while cooling is shown in Figure [2j The data are pre- 
sented as p(t)/p(t = 0) so that the values are normalized 
to unity at t = for easy comparison. We found that 
that below 160 K, the resistivity of the sample increases 
with time considerably. A maximum relative increase 
in resistivity of about 200 % for a duration of one hour 
is seen at 145 K, the time dependence being lower both 
above and below this temperature. We fitted the p(T, t) 
curves in figure 2 to the stretched exponential function 

p(t) = p + Pl {l-e-(iy) (2) 

where po, pi, r and 7 are fit parameters. The fits 
are quite good with the R 2 value greater than 0.999 in 
most cases. See Table [H We note that the exponent 7 
lies in the range 0.5 < 7 < 0.6 and r has a peak around 
147.5 K. The variation of resistivity with time shows that 
the system slowly evolves towards an insulating state at a 
constant temperature. We collected data up to 12 hours 
(not shown here) to check whether the system reaches 
an equilibrium state, but found that it was continuing to 
relax even after such a long time. 

The magnitude of time dependence in heating runs 
(maximum m 0.2%) is negligible compared to what one 
gets in cooling runs (maximums 200%) (figure [2]), which 
suggest that in the heating run, below Tmi, the sample 
is almost fully insulating and stable. 

Figure [3] compares the resistivity in a cooling run with 
and without intermediate aging. These data were taken 
as follows: we start from 220 K, come down to 160 K at 
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1.22(1) 


4.04(6) 


0.582(2) 


0.0014 


0.99988 
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0.391(5) 


2.6(1) 


0.557(7) 


0.0003 


0.99845 



Table I: Fit parameters for the time dependence data shown 
in Figure [2] The degrees of freedom of the fits DOF w 1000. 
The x 2 /DOF for 150 K and 155 K are too small, indicating 
that we have overestimated the error in resistivity in these 
cases. Anyway, we note that, the R? values are consistently 
good and indicate reasonably good fits. 
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Figure 3: (Color online) Temperature dependence of p on 
cooling with ageing by intermediate stops of one hour each 
(blue open circles) at 155, 150, 145, 140, 135, and 130 K. The 
black line show the resistivity without ageing. The red line 
shows the results of the simulation discussed in section IIVDI 



2K/min, collect time dependence data for one hour and 
after that resume cooling at 2K/min and go down to 
155 K, collect resistivity time dependence again for one 
hour and so on down to 110 K for every 5K interval. 
We note that when cooling is resumed after aging for 
an hour, p(T) curve merges smoothly with the curve ob- 
tained without aging within a few kelvin. A rather simi- 
lar observation has been reported in the phase-separated 
manganite Lao.5Cao.5Mno.95Feo.05O3.- 



IV. DISCUSSION 

We have seen that, below T MI , while cooling, the sys- 
tem is not in equilibrium and evolves with time which 
suggests that it is in a metastable state. Or, possibly, it 



consists of at least one metastable constituent. We infer 
that the resistivity of this metastable constituent slowly 
increases with time which probably implies that a slow 
metal to insulator transition is going on in the system. 
On the other hand while heating up from low tempera- 
ture we saw that the system remains insulating all the 
way up to 188 K and then it transforms to the metallic 
state by 200 K. The resistivity was found to have negli- 
gible time dependence during the whole of the heating 
run. This result suggests that while heating from low 
temperature the system is in a stable equilibrium state. 

It is well known that below a phase transition tem- 
perature a high temperature phase can survive as a 
metastable supercooled state. Based on this information 
we propose that below the M-I transition temperature, in 
the hysteresis region, our system consists of supercooled 
(SC) metallic regions and stable insulating regions. A 
supercooled metallic region would be separated from its 
stable insulating state by an energy barrier which can 
be crossed with the help of perturbations such as ther- 
mal fluctuations, mechanical disturbances and so on. The 
crossing of the barrier from the metallic state into a stable 
insulating state can give rise to time dependence of resis- 
tivity in the system. Now, it is known that, on cooling a 
supercooled system it can either switch to a stable state 
or get kinetically arrested into a glassy state. From the 
fact that no time dependence is observed while heating 
from low temperature towards the M-I transition tem- 
perature we rule out the possibility of any glassy phases 
forming in the system. For, if a glassy phase did form, 
time dependence would have been observed while heat- 
ing from low temperature because of the de-arrest of the 
same. Consequently we are forced to conclude that the 
supercooled phases which are present while cooling would 
switch to the stable insulating state at a sufficiently low 
temperature. At a low enough temperature all of the ma- 
terial would be in the stable insulating state and hence 
subsequent heating will not show any time dependence. 

Based on the fact that no time dependence is observed 
above the M-I transition temperature, while heating, we 
rule out the possibility of the existence of superheated 
phases in our system. 

In the next subsections we get down to the nitty-gritty 
of the model we use to understand the experimental ob- 
servations. 



A. The Model 

In a first order transition, a metastable SC phase can 
survive below the first order transition temperature (Tc), 
till a certain temperature called the limit of metasta- 
bility (T*) is reached . 3 i 31 i 32 In the temperature range 
T* < T < Tc there is an energy barrier separating the 
SC metastable phase from the stable insulating state. 
The height of the energy barrier, U, can be written as 
U oc f(T — T*), where / is a continuous and monotonic 
function of (T — T*), and vanishes for T < T*. As the 
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temperature is lowered, at T — T*, the SC metastable 
phase becomes unstable and switches over to the stable 
insulating state as the energy barrier becomes zero in this 
case^i At T > T* the SC metastable phase can cross over 
to the stable insulating state with a probability (p) which 
is governed by the Arrhenius equation 

p oc e" 7 ^ (3) 

which tells us that the barrier will be crossed with an 
ensemble average time constant r oc 1/p. If we imagine 
an ensemble of such SC phases with the same barrier U, 
then the volume of the metastable phase will exponen- 
tially decay with a time constant r. 

A transition from a supercooled state to a stable state 
is an avalanche-like transition which happens abruptly. 
Thus a single crystal or a crystallite, in the case of a 
polycrystalline material, would remain in the SC state 
above its temperature of metastability, and it will switch 
to the stable state as a whole when it is pushed over the 
free energy barrier by an energy fluctuation or if in the 
process of cooling the temperature of metastability is at- 
tained. Examples for this would be the sharp M-I transi- 
tions along with hysteresis and time dependence observed 
in VO, V 2 3 , Li and Na&±£>±& It has been claimed that 
the presence of defects, strains and non-stoichiometry can 
give rise to different T*'s for different crystallites, or it 
may even produce regions having different T*'s inside a 
single crystal or crystallite . 3 ' 33 ' 34 ' 35 We will call a crys- 
tallite or a region within a crystallite that will switch as a 
singe entity as a switchable region (SR) in our discussion. 

Carrying forward the above arguments to our poly- 
crystalline system, which is made up of tiny crystallites, 
we can say that it will be made up of a large collec- 
tion of SR's. Each SR will have a unique temperature of 
metastability, T*, and a volume, V. The energy barrier 
U that we discussed earlier, being an extensive quantity, 
will be proportional to the volume of the SR and we can 
write 

U = Vf(T - T*) (4) 

where / is the continuous and monotonic function which 
vanishes for non-positive values of its argument. This 
means that the various SR's with their different T* and 
V will have different energy barriers, which implies that 
the time constant r will spread out and become a distri- 
bution of time constants depending on the distribution 
of the size and T* of the SR's. This can give rise to 
the volume of the metallic state decaying in a stretched 
exponential manner with time i 36 ' 37 This behavior of the 
metallic volume with time will lead to the resistivity also 
evolving with time in a similar fashion. We shall see in 
the next subsection how the metallic volume and the re- 
sistivity are related to each other. 
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Figure 4: (Color online)Temperature variation of V (= 1 — /), 
the insulating volume fraction. The blue circles show the cool- 
ing cycle (2 K/min cooling rate) and the red squares represent 
the heating cycle. The green pluses represent the insulating 
volume fraction for 0.2 K/min cooling rate. The dashed hori- 
zontal line with the label Vth(= 84% insulating volume) rep- 
resents the percolation threshold. The inset shows the vari- 
ation of \dV/dT\ for the cooling cycle (blue circles: cooling 
rate 2 K/min, and green pluses: cooling rate 0.2 K/min). 

B. Volume Fraction Calculation 

To check the validity of our model through numerical 
simulations we need to deal with the volumes of super- 
cooled and stable phases. Our experimental measure- 
ments are of resistivity and we need to find a way of esti- 
mating the volumes of the different phases from our data. 
It is well known that the conductivity of a binary mix- 
ture consisting of insulating and metallic phases depends 
on their respective volume fractions, geometries, distri- 
bution, and conductivities of the insulating and metallic 
phases (07 and cm)- 

We shall use an easy to use formula given by 
McLachlan~ which is based on a general effective 
medium (GEM) theory, for doing this conversion. The 
McLachlan GEM equation has been successfully applied 
to a wide variety of isotropic, binary, macroscopic mix- 
tures and it has been seen to work well even close to the 
percolation threshold i 38 ' 39 ' 40 If <je is the effective electri- 
cal conductivity of a binary MI mixture, the GEM equa- 
tion says 

1 V* l/t\ 1 lit l/u 

a f] K - zA ) , f « z zA ) - ^ 

(a/ +Aa^ ) (cr^ +AaJ ) 
where / is the volume fraction of the metallic phases and 
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A = (1 — f c )/fc, fc being the volume fraction of metallic 
phases at the percolation threshold, and t is a critical 
exponent which is close to 2 in three dimensions . 39 ' 41 
The constant f c depends on the lattice dimensionality, 
and for 3D its value is 0.164^ 

In order to calculate the volume fraction of metallic 
and insulating phases from Equation |(5]), we need their 
respective resistivities pu and pi as functions of temper- 
ature. Pm{T) was obtained using pu = Po + /3T, where 
/? is the temperature coefficient of resistivity estimated 
from the resistivity data above the M-I transition. pi(T) 
was calculated using Equation (U) with the parameters 
obtained by fitting the resistivity data below 115 K. Us- 
ing the above information and the resistivity data of the 
major loop shown in Figure [IJ we calculated the volume 
fraction of metallic and insulating phases and it is shown 
in Figure |H In the cooling cycle the volume fraction of 
the insulating phases (V = 1 — /) slowly increases on de- 
creasing the temperature below Tmi , while in the heating 
cycle it remains nearly constant up to about 185 K, and 
then drops to zero by 200 K. As can be seen from Figure 
|4]the percolation threshold for the cooling runs occurs at 
around 144 K. Below this temperature there will be no 
continuous metallic paths in the system. 
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Figure 5: (Color online) AV vs time for 140 K to 155 K. The 
gray symbols show the values derived from experiments and 
the black lines show the simulated values. The inset shows 
how the parameter a depends on temperature. The red lines 
indicate the values we have chosen for the simulation. 



Figure [5] displays the increment of insulating volume 
fraction, AV(exp), as a function of time that has been 
extracted from the time dependent resistivity data of Fig- 
ure [2 We see that the maximum time dependence in 
volume fraction is seen at 147.5 K which is not the same 
temperature at which the maximum time dependence in 
resistivity is seen. It is interesting to note that the rate of 
change of metallic volume fraction, \dV/dT\, has a maxi- 
mum around 147.5 K in the cooling runs (inset of Figure 
H]) and it coincides with the temperature of the maximum 
volume time dependence. 

During cooling \dV/dT\ of the inset of Figure H] rep- 
resents the amount of volume that will change from the 
supercooled metallic to the insulating state for a unit 
temperature change. Now a small change in temperature 
from T to T — ST will result in all the supercooled SR's 
with their T* falling in that temperature range switching 
to the insulating state from the supercooled state. This 
means that \dV/dT\8T at T is a good measure of the 
volume fraction of supercooled metastable regions which 
have their T* close to T within a temperature range ST. 
Thus \dV/dT\ of the cooling curve represents the volume 
distribution of T*'s in the system. In coming to the above 
conclusion we have disregarded the small fraction of SR's 
that would be switching due to the time elapsed in cov- 
ering the small temperature change. The justification for 
this is the fact that the \dV/dT\ values calculated from 
the 2 K/min and 0.2 K/min cooling curves are practically 
indistinguishable as can be seen from the inset of Figure 
H 



C. Details of the Simulation 

We carried out Monte Carlo simulations to try to un- 
derstand the experimental data. We take a distribution 
of 10 5 SR's which have their volumes (V) uniformly dis- 
tributed in a certain range (Vo,Vmax)- These SR's were 
assigned a T* in such a way that the volume distribu- 
tion of the supercooled metallic SR's matches with the 
\dV/dT\ curve shown in the inset of Figure HI These 
metallic SR's have an energy barrier given by Equation 
(j4]) that separates them from the insulating equilibrium 
state. Now we assume that f(T — T*) has the simple 
power law form / = c(T - T*) a for T > T*, where c and 
a are positive constants, and / = for T < T* . Using 
Equation O at any given temperature, the probability of 
switching from the metallic to the insulating state in a 
certain time interval is taken to be 

p = e - aV — t — (6) 

with a = c/kB- 43 

Now to carry forward the simulation we have to esti- 
mate the parameters Vq, V max , a, and a. We decided to 
first simulate the time dependence of the insulating vol- 
ume evolution shown in Figure [5] for this purpose. The 
time dependence of the insulating volume is generated 
as follows. To begin with we quench the simulated sam- 
ple from above Tmi to the temperature of interest. In 
this process all the SR's which have their T*'s above the 
temperature of interest are switched from the metallic to 
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the insulating state. Now the switching probability of 
each of the remaining metallic SR's is compared with a 
uniform random number in the interval [0,1) and if the 
switching probability is greater than the random num- 
ber, the metallic SR is flipped to the insulating state and 
the increase in the insulating volume fraction is recorded. 
This operation is carried out for each one of the metal- 
lic SR's. The time required for this process is taken to 
be the time interval at which the experimental data was 
recorded, which in our case is 3.6 seconds. We repeat 
this process one thousand (n) times to generate data for 
one hour which is the time over which the experimental 
data was collected. The parameter space of Vq, V max , a, 
and a was searched to minimize the error function 



Err(Vo,V max ,a, a) = I (&V{exp) l - AV(cal)i) 

1=0 

(7) 

where AV(exp)i is the increase in insulating volume frac- 
tion obtained from experimental data, AV(cal)i is the in- 
crease in insulating volume fraction calculated from the 
simulation and the subscript i stands for the number of 
the time step. We found that as far as the volume of the 
SR's are concerned it is the ratio of V max to Vb that de- 
cides the minimum value of Err(Vo,V ma x,a,a) and not 
the individual values. So we fixed Vb to be equal to unity 
and found that overall V max ~ 6 and a « 0.25 gave the 
best values for the error function with parameter a lying 
in the range 270 to 400 depending on the temperature. 
We refined the simulation by fixing the values of V max to 
be 6 and a to be 0.25 and varying only a. We compare 
in Figure [5] the experimental data (gray symbols) and 
the simulation (black lines) and see that the simulation 
reproduces the experimental data quite closely. In the 
inset of the figure we show how the parameter a varies 
with temperature.— For further calculations we fixed the 
value of a to be 280 up to 146 K and above that temper- 
ature a was taken to vary linearly with temperature as 
a = 280 + 14(T — 146) and this is shown as red lines in 
the inset of Figure [H 

It would have been ideal if all the parameters of the 
simulation turned out to be constants. But we find that 
we have to introduce some temperature dependence in 
a to simulate the experimental results accurately. The 
temperature dependence of a could possibly be related 
to the percolation threshold which is at 144 K. Below 
this temperature the background matrix for most of the 
SR's would be insulating while above this temperature 
the background would gradually change to a metallic one 
as one goes away from the percolation threshold. We 
also note that an SR in its insulating state would have 
a slightly more distorted crystal structure along with a 
somewhat larger volume. Thus we see that the surround- 
ings of an SR would be very different above and below 
144 K and this may have a bearing on the energy barrier 
the SR has to cross and this could result in a temperature 
dependent a. 



In the simulation we quench the sample from above 
Tmi to the temperature of interest while in the ex- 
periment the sample is cooled at a fixed cooling rate 
(2K/min) to the temperature of interest. In the sim- 
ulation we can cool the sample slowly only if we know 
the parameters, and since we have no idea of the param- 
eters to begin with, as a way out we decided to quench 
the sample. As we shall see soon, the simulation and 
the experimental data agree quite well for the other ex- 
periments and hence it was felt that further refining of 
the simulation parameters would not improve the simu- 
lation very much and hence we decided to stick with the 
parameters obtained by quenching. 

D. Simulation Results 

We will now describe the simulation results on resistiv- 
ity hysteresis and minor loops, cooling rate dependence 
and rejuvenation. 



1. Hysteresis and Minor Loops in Resistivity 

Hysteresis and minor loops in resistivity are calculated 
by first simulating how the insulating volume evolves as 
one cools and heats the (simulated) sample. The cool- 
ing was done by lowering the temperature by 0.12 K first, 
during which process all those SR's with their T* in that 
0.12 K range are switched to the insulating state, fol- 
lowed by a time step of 3.6 seconds, during which all 
the remaining metallic SR's are given a chance to relax 
to the insulating state. We consider the above proce- 
dure to be a good approximation for cooling at a rate of 
2K/min. This process is repeated as many times as re- 
quired to complete the cooling part of the simulation. In 
the heating part of the simulation we again take temper- 
ature steps and time steps as before. But it may be noted 
that a temperature step while heating will not switch any 
of the SR's to the insulating state because all those SR's 
with T*'s falling inside the temperature step would have 
already been switched to the insulating state during the 
preceding cooling process. During a time step all of the 
metallic SR's do get a chance to switch their states, but 
we point out that during heating the chance of a metallic 
SR switching to the insulating state would decrease dras- 
tically with increasing temperature because the energy 
barrier is a monotonically increasing function of T — T* 
and all the remaining metallic SR's would have their T* 
below the turning point, i.e., the temperature at which 
the cooling was stopped and heating started. After the 
volume simulation is completed we translate the results 
back to resistivity using Equation © . 

In Figure[6]we show the simulated results for hysteresis 
and minor loops and we note that they are very similar to 
the experimental results shown in FigureQ] except for the 
sharp transition in the simulation during heating. This 
can be attributed to the fact that we have not taken into 
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Figure 6: (Color online) Simulated results for the hysteresis 
along with minor loops. The labeling and color coding are the 
same as in Figured] Cooling rate dependence of resistivity is 
shown in the inset where the symbols represent the experi- 
mental data and the lines are from simulation. Blue star and 
blue line: 2K/min, red circles and red line: 0.2K/min, black 
line: 1 K/hour simulation, black pluses: infinitely slow cooling 
rate, estimated from experiment. 



account the broadening of the M-I transition due to dis- 
order and finite size effects in our simulation. It is inter- 
esting to note that the simulated minor loops reproduce 
quite closely some of the features of the experimental 
loops such as the positive slope of the heating (top) part 
of loop a and the negative slope at the beginning of the 
heating part of loop c. The heating part of curve a has 
a positive slope because in this case most of the SR's are 
in their insulating state and the temperature dependence 
of their resistivity is the dominating factor which deter- 
mines its behavior. On the other hand the beginning of 
the heating part of curve c has a negative slope because 
in this case most of SR's are in the metastable metallic 
state and the increase of their resistivity with time is the 
dominating factor here. 



2. Cooling Rate Dependence 

Simulating the cooling rate dependence is straightfor- 
ward. The time steps were fixed at 3.6 seconds and the 
temperature steps were adjusted in size to simulate faster 
or slower cooling. The results are shown in the inset of 
Figure[6l It is clear that the experimental data and the 
simulation agree very well for 0.2K/min and 2K/min 
cooling rates. We have also shown the simulation result 
for IK /hour cooling rate and it is seen to fall in between 



the 2 K /man data and the infinitely slow data. It is easy 
to see that a slower cooling rate will result in a higher 
resistivity because more of the SR's switch to the insu- 
lating state during slow cooling thus resulting in a higher 
resistivity. 



3. Intermediate Ageing 

Intermediate ageing is also fairly easy to simulate. We 
cool the simulated sample at 2 K/min to the temperature 
of interest, do one thousand time steps to simulate wait- 
ing for one hour, and resume cooling and so on. In the 
inset of Figure [3] we have shown the simulation results 
for intermediate ageing and we see that the agreement 
between data and simulation is reasonable. 

From Figure[3]it is clear that on resuming cooling after 
an intermediate stop of one hour the cooling curve merges 
with the curve obtained without ageing within about 3 K 
or less. When we stop the cooling and age the sample at 
a fixed temperature the supercooled metallic SR's with 
relatively small U will switch over to the insulating state. 
As can be inferred from Equation Q those SR's will have 
a small U which have (i) their metastability temperature 
(T*) close to their temperature (T) (ii) a small size. The 
major contributions to resistivity change will come from 
the relatively larger SR's making the transition from the 
metallic state to the insulating state. In the light of this 
the merger of the cooling curves with and without ageing 
within a small temperature change of about 3 K means 
that most of the larger SR's which undergo the transi- 
tion from metal to insulator during ageing have their T* 
within a few kelvin of the temperature of the sample. 



E. Predictive Power of the Model 

A model is no good if it cannot predict new results. So 
to test our model we decided to do a simulation first and 
then see whether experiment will reproduce it. In the 
simulation we oscillate the temperature of the sample 
between 145 K and 146 K for one hour immediately after 
cooling it from 220 K to 145 K. All temperature variations 
were done at 2 K/min. Figure[7fa) shows the results of 
the simulation. We see that the resistivity keeps on in- 
creasing on repeated thermal cycling; the increase per cy- 
cle is larger at the beginning and slowly it tapers off. We 
repeated the simulation on the heating cycle also, cooling 
the sample first from 220 K down to 85 K and then bring- 
ing it back up to 145 K before starting the temperature 
oscillations. No detectable change in resistivity was seen 
with thermal cycling. 

We expected that during the cooling cycle the temper- 
ature oscillations would disturb the sample and we would 
get a larger increase in resistivity than if the temperature 
were kept fixed at 145 K. See Figured But instead, we 
find intriguingly, that the increase in resistivity while os- 
cillating the temperature is considerably less than what 
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Figure 7: (a) shows the simulation results of how the resistiv- 
ity evolves if the temperature of the sample is oscillated be- 
tween 145 K and 146 K for one hour after cooling from above 
Tmi- Po is the resistivity of the sample when it reaches 145 K 
after cooling from 220 K at 2K/min. The arrows indicate 
that the resistivity keeps on increasing on thermal cycling, 
(b) shows the experimental results for the same. 

one gets at the constant temperature of 145 K. Even more 
surprisingly it is less than that observed at 147.5 K. It will 
be very interesting to see if experiment agrees with this 
prediction. 

In Figure [7]Jb) we show the experimental results ob- 
tained during thermal cycling after cooling the sample 
from 205 K. One can see that the agreement between the 
simulation and the experiment is very good. We repeated 
the thermal cycling after cooling the sample down to 85 K 
and bringing the temperature back up to 145 K, just as 
in the simulation. Here also the experiment agreed very 
well with the simulation; no detectable change in resis- 
tivity was found. 

The unexpectedly lower resistance rise when the tem- 
perature is oscillated after cooling can be attributed to 
(1) the smaller number of SR's available for switching 
because all those SR's with T* falling in the range 145 K 
to 146 K would have already switched to the stable insu- 
lating state the first time the temperature was lowered 



to 145 K and (2) a higher average barrier height which 
an SR has to overcome, because the barrier height goes 
as(T-T*) a . 



V. CONCLUSION 

Our experimental results and its excellent agreement 
with the simulation suggests strongly that while cool- 
ing the physical state of NdNiOa is phase separated be- 
low the M-I transition temperature; the phase separated 
state consists of SC metallic and insulating regions. A 
metastable metallic region switches from the metallic to 
the insulating state stochastically depending on the close- 
ness of the temperature of metastability and the size of 
the region. At low temperature, below 115 K or so, the 
system is insulating, all the SR's having switched over 
to the insulating state. While heating the SR's remain 
in the stable insulating state till the M-I transition tem- 
perature is reached and then switch over to the metallic 
state. 

In our study we have developed a physical explanation 
of the time dependence effects observed in NdNiOa based 
on phase separation and supercooling. We believe that 
the physics developed here may have implications for first 
order solid-solid phase transitions in general. We have 
shown that our model has predictive power and we feel 
that it may serve as a useful template for understanding 
slow dynamics in hard condensed matter in such cases 
where there is phase separation. The apparently counter 
intuitive simulation prediction in the temperature oscilla- 
tion experiment and its experimental confirmation gives 
us confidence that our model has captured the underlying 
physics of the problem truthfully. 
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